USITP/98-08 
MPI-PhT/98-43 
June 1998 



Clumpy Neutralino Dark Matter 



Lars Bergstrom 

Department of Physics, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden; 

lbe@physto.se 



Joakim Edsjo 

Center for Particle Astrophysics, University of California, 
301 Le Conte Hall, Berkeley, CA 94720-7304, U.S. A; 
edsjo@cfpa. berkeley. edu 



Paolo Gondolo 

Max Planck Institut fiir Physik, Fdhringer Ring 6, 80805 Munich, Cermany; 

gondolo @mppmu. mpg. de 



Piero Ullio 

Department of Physics, Stockholm University, Box 6730, SE-113 85 Stockholm, Sweden; 

piero@physto.se 



Abstract 



We investigate the possibility to detect neutralino dark matter in a sce- 
nario in which the galactic dark halo is clumpy. We find that under customary 
assumptions on various astrophysical parameters, the antiproton and contin- 
uum 7-ray signals from neutralino annihilation in the halo put the strongest 
limits on the dumpiness of a neutralino halo. We argue that indirect detection 
through neutrinos from the Earth and the Sun should not be much affected 
by dumpiness. We identify situations in parameter space where the 7-ray 
line, positron and diffuse neutrino signals from annihilations in the halo may 
provide interesting signals in upcoming detectors. 







I. INTRODUCTION 



The mystery of the dark matter in the Universe remains unsolved. Among the more 
plausible candidates (not only needed to solve the dark matter problem) can be found the 
neutralino, the lightest supersymmetric particle in the Minimal Supersymmetric Standard 
Model (MSSM) (for a review, see Another candidate is for instance the axion which is 
still a viable option for a narrow range of axion masses 0] . Irrespective of the exact nature 
of the dark matter, there are reasons to believe that its distribution in the dark halos of 
galaxies need not be perfectly smooth For instance, early fluctuations in the dark 

matter may go non-linear long before photon decoupling, evading the argument of slow, 
linear growth after recombination. Also, if cosmic strings or other defects exist, they may 
seed the formation of density-enhanced dark matter clumps. 

Since very little is known about the inherently non-linear problem of generating dark 
matter clumps, in this paper we will use a phenomenological approach where we simply 
assume the existence of clumps with a given density profile, making up a certain fraction 
of the total mass of the Milky Way halo. We investigate the effect of this dumpiness on 
the various proposed detection methods for neutralino dark matter. Detection rates depend 
crucially on the neutralino distribution in momentum and position space. Some detection 
rates, in particular those of antiprotons and photons generated by neutralino annihilations 
in the galactic halo, increase substantially compared to the case of a smooth dark matter 
distribution. For a given set of parameters of the supersymmetric models (such as mass and 
couplings of the neutralinos) we can then use present experimental limits on these fluxes to 
bound the degree of dumpiness allowed in that particular dark matter model. Alternatively, 
given a positive experimental signature, we can identify regions in the combined parameter 
space of halo dark matter distribution and supersymmetric models to identify candidates 
consistent with the data. This approach was used recently by three of us in connection with 
new data from the EGRET gamma ray detector . Some of our results may be of interest 
also in the standard non-clumpy scenario, which is of course included in our treatment and 
is easily recovered by putting the fraction of the halo in the form of clumps equal to zero. 

II. CLUMPINESS IN THE MILKY WAY HALO 

Present observational data give very poor constraints on the distribution of dark matter 
in the galaxy. The dynamics of the outer satellites of the galaxy clearly indicates that 
luminous matter provides just a fraction of the total mass of the Milky Way and that the 
major contribution must come from a dark matter halo whose size is larger than the radius 
of the disk. Nevertheless it is not possible to extract from present kinematic information any 
accurate knowledge of the density profile of the dark matter halo. It is, however, natural 
to assume that galactic dark matter profiles obey a law of universality. Then, a possible 
approach is to infer the functional form of the Milky Way halo density profile from the 
results of A^-body simulations of hierarchical clustering in cold dark matter cosmologies, 
fitting the normalization parameters to known dynamical constraints. This approach has 
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been followed in Ref . M : among the general family of spherical density profiles, 
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(1) 



it was considered the case of the Kravtsov et al. profile ^ which is mildly singular towards 
the galactic centre with 7 ~ 0.2-0.4, of the Navarro et al. profile which is more cuspy 
(7 = 1), and, for comparison the modified isothermal distribution, 7) = (2,2,0), 

extensively used in dark matter detection computations. 

The dark matter density profile inferred in this way should be regarded as the function 
that describes the average distribution of dark matter in the galactic halo; the standard 
assumption which is generally made at this stage is that dark matter particles in the halo 
form a perfectly smooth 'gas'. This approach is in some way arbitrary: although the dark 
matter particle distribution has to be regarded as smooth on intermediate length scales, 
probably around 0.01-1 kpc, there are reasons to question whether this is true on smaller 
scales. We here entertain the possibility that at least a fraction of the dark matter in the 
halo is clustered in substructures with high matter density, 'clumps' of dark matter. Several 
authors have introduced dumpiness as a generic feature of cold dark matter cosmologies. Silk 
and Stebbins |Q have considered clump formation in cosmic string, texture and inflationary 
models, giving also predictions for survival to tidal disruption (see also Ref. [|]). Kolb 
and Tkachev |^ have studied isothermal fluctuations giving very high-density dark matter 
clumps. 

Simulations of structure formation in the early Universe do not yet have the dynamical 
range to give predictions for the size and density distribution of small mass clumps (we 
focus here mainly on clumps of less than around 10^ solar masses which avoid the problem 
of unacceptably heating the disk The formation of clumps on all scales is however a 

generic feature of cold dark matter models which have power on all length scales. If self- 
similarity is a guide, galaxy halos may form hierarchically in a similar way to that of cluster 
halos (see e.g. Ref. [p!0| ]). 

Rather than examining the different scenarios for clump formation, we take a more phe- 
nomenological approach and perform a detailed discussion on the implications of dumpiness 
on neutralino dark matter searches. 

We thus simply postulate that a fraction / of the total dark matter is concentrated in 
clumps, which are assumed to be spherical bodies of typical mass Md and matter density 
profile Pci{f\:i)- The total number of clumps inside the halo is given by: 

iV. ~ ^ (2) 

where is the total mass of the halo. Two opposite scenarios seem to be plausible. There 
might be few heavy clumps, with masses up to maybe M^i ~ 10^-10*Mq (above which 
the local gravitational distortion effects would be too severe); such massive bodies could 
in principle be identified from the analysis of the rotation curves of the galaxy, they may 
however have escaped observation so far because their detection in this way may be difficult 
if the fraction of dark matter in clumps is small, say / ~ 1%. A second possibility is that 
clumps are much lighter, with M^i less than lO^-lO^M©, in which case larger fractions of the 
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TABLE I. The ranges of parameter values used in our scans of the MSSM parameter space. 
Note that several special scans aimed at interesting regions of the parameter space have been 
performed. In total we have generated about 85000 models that are not excluded by accelerator 
searches. 



halo mass might be in clumps (in the extreme scenario all of it) . In the many small clumps 
scenario, on which we mainly focus, we can define a probability density distribution of the 
clumps in the galaxy which in the limit of large /, to fulfill dynamical constraints, has to 
follow the mass distribution in the halo. 

Consider a Cartesian coordinate system with origin in the galactic centre. Then the 
probability for a given clump being in the volume element d^x at position x is 

Pci{x) d^x = ^ p{x) d^x (3) 

which has the correct normalization / Pci{x)d^x = 1. 

It is convenient to introduce the dimensionless parameter 6 

r 1 Jd^r^i (Pc/(rd))^ , . 

= ^ 4 

Po J d-^Tci Pci{rci) 

which gives the effective contrast between the dark matter density in clumps and the local 
halo density po- For a dark matter density inside the clumps which is roughly constant, pd, 
it reduces to the form 

5 = ^. (5) 

Po 

We show in Section that in the many-clumps scenario it is just the product f6 which 
determines the increase of the signal compared to a smooth halo in most indirect detection 
methods. The product f6 is directly related to the ratio of the total dark mass in clumps 
to the volume of a typical clump. 



III. SUPERSYMMETRIC MODELS FOR DARK MATTER 

Although it is possible that the halo dark matter may be composed of particles not yet 
predicted by particle physics models, it is very attractive to assume that they are weakly 
interacting massive particles (WIMPs). Massive particles with weak interactions give a 
relic density which is of the right order of magnitude to explain the dark matter on all 
scales from dwarf galaxies and upwards. We will consider a specific class of such particles, 
supersymmetric particles, which is general enough to illustrate the effects of dumpiness. 
Our results should be of more general validity, however. 
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We work in the Minimal Supersymmetric Standard Model (MSSM) as defined in 
Refs. [|ri| l|. For details on our notation, see Ref. [12|. The lightest stable supersymmetric 



particle is in most models the neutralino, which is a superposition of the superpartners of 
the gauge and Higgs fields, 

It is convenient to define the gaugino fraction of the lightest neutralino, 

= |iVn|' + |iVi2|' (7) 
For the masses of the neutralinos and charginos we use the one-loop corrections as given in 



The MSSM has many free parameters, but with some simplifying assumptions, we are left 
with 7 parameters, which we vary between generous bounds. The ranges for the parameters 
are shown in Table |. For the detection rates of neutralino dark matter we have used the 
rates as calculated in Refs. ||6|J^,p!^-[T7 . 



We will throughout this paper assume that the neutralinos make up most of the dark 
matter in our galaxy. We only consider therefore MSSM models which are cosmologically 
interesting, i.e. where the neutralinos can make up a major fraction of the dark matter in 
the Universe without overdosing it. We will choose this range to be 0.025 < VL^Ji^ < 0.5. 
For the relic density calculations we have used the detailed calculations performed in Ref. 

0- 



IV. DETECTION METHODS CONSTRAINING CLUMPINESS 

Some observational consequences of a clumpy dark matter halo have been pointed out 
previously, such as the obvious gain in gamma ray signal from annihilation in the halo since 
the flux from a particular volume element is proportional to the square of the dark matter 
density there P,p|,p!^-plll • Also, in Ref. it was noted that the antiproton flux could be 



enhanced, although the treatment was sketchy and not entirely correct concerning the way 



the rescaling was done. In Ref. it was investigated whether encounters with dark matter 
clumps on geophysical time scales could have left imprints in ancient mica. 

As we show in this section, indirect detection through cosmic antiprotons and gamma 
rays set the most stringent limits on clumpy neutralino dark matter, therefore we investigate 
these cases flrst. 



A. Gamma-rays 

Since gamma rays produced in neutralino annihilations in the halo travel in straight 
paths essentially without any absorption, and since the annihilation rate and hence the flux 
would be enhanced by clumps along a particular line-of-sight, the effects of dumpiness are 
easy to understand. 

Neutralino annihilation in the galactic halo may produce both a 7-ray flux with a con- 
tinuum energy spectrum and monochromatic 7-ray lines. 
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The continuum contribution (see Ref. and references therein) is mainly due to the 
decay of 7r° mesons produced in jets from neutrahno annihilations. To model the fragmenta- 
tion process and extract information on the number and energy spectrum of the 7s produced 



we have used the Lund Monte Carlo Pythia 6.115 |2J]. We have performed the simulation 
for 18 neutralino masses between 10 and 5000 GeV and for the cc, bb, ti, W^W^, Z^Z^ 
and gg annihilation states. For each final state and for each neutralino mass we have simu- 
lated 2.5 X 10^ events which are tabulated logarithmically in energy. For any given MSSM 
model, we then sum over the annihilation channels and interpolate in these tables. For 
the annihilation channels not included in the simulations, like the ones with one gauge and 
one Higgs boson as well as those with two Higgs bosons the flux is calculated in terms of 
the flux from the simulated channels. We include all two-body final states at the tree level 
(except light quarks and leptons) and the one-loop processes Zj and gg. For final states 
with Higgs bosons, we let the Higgs bosons decay in flight by summing the contributions to 
the gamma flux from the Higgs decay products in the Higgs rest system and then boost the 
spectrum averaging over decay angles. Given the annihilation branching ratios we then get 
the spectrum for any given MSSM model. The continuum signal lacks distinctive features 
and it might be difficult to discriminate from other possible sources. It will however be a 
powerful tool to put constraints on the dumpiness parameters. 

A much better signature than the continuum contribution is given by monochromatic 
7-ray lines which arise from the loop-induced S-wave neutralino annihilations into the 27 
and Z7 final states and which have no conceivable background from known astrophysical 
sources. The amplitude of these two processes in the MSSM was computed only recently at 
full one loop level p5| , p6| . Large deviations from previous partial results (see Ref. |]1| and 
references therein) were found, in particular it was pointed out that a pure heavy Higgsino 
has a remarkably high annihilation branching ratio both into 27 and Z'j, adding at least 
a factor of 10 to previous estimates of the 27 line. A detailed phenomenological study is 
given in Ref. 0] where a smooth halo scenario was considered and it was shown that the 
monochromatic lines could be detected by the new generation of space- and ground-based 7- 
ray experiments, provided that a sensible enhancement of the dark matter density is present 
towards the galactic centre. We examine here the perspectives of detecting the continuum 
and the line signals in a given clumpy scenario. 

Consider a detector with an angular acceptance AQ pointing in a direction of galactic 
longitude and latitude {i,b). The gamma ray fiux from neutralino annihilations at a given 
energy E is given by 

$^(E, Afi, i, b) ^ 1.87 ■ 10-s — ( J (£, 6) ) (A(]) cm-2 s'^ sr'^ . (8) 

dE 

In this formula we have defined a factor dS/dE which depends on the nature of relic WIMPs. 
For the continuum 7-ray signal, the 27 line and the Z7 line signal, respectively, it is given 
by: 




10GeV\ ^/ vap \ dN, 



E 



J ^ ^10"^^ cm3 s-V dE 



F 



7 



lOGeVX / 2w2, 



^7 



J Vl0-26cm3s-i 



5 



Here is the neutralino mass, F are the allowed final states which contribute to the 
continuum signal as specified above. For each of these, vap is the annihilation rate and 
dN^ /dE is the differential energy distribution of produced photons. The product of relative 
velocity and cross section va2'y is the annihilation rate into the 27 final state (as given in 
Ref. p5[). Similarly, vaz-y is the rate into the Z7 final state (as given in Ref. ||26|). In Eq. (H) 



the dependence of the fiux on the dark matter distribution, the direction of observation (£, h) 
and the angular acceptance of the detector Af2 is contained in the factor ( J (£, h) ) ( AQ) . 
If we assume a spherical dark matter halo in the form of a perfectly smooth distribution of 
neutralinos, it is equal to 

( J (V^) ) (Afi) = f dn' [ dL ( — . (10) 

^ ^^^^^ ' 8.5kpcAfiiAn June of sight VO-3 GeV/cm^/ ^ ^ 

Here L is the distance from the detector along the line of sight, ip is the angle between the 
direction of observation and that of the galactic center, related to {i, h) by cos tp = cos i cos b. 
The integration over dQ' is performed over the solid angle AQ centered on ip. 

We shall now examine the consequences of introducing clumps in the halo. The contin- 
uum 7-ray signal in the few clumps scenario has as mentioned been examined in some detail 
in the literature. Following the approach of Refs. |ll8| , [l9| and estimating the most probable 



distance for the nearest clump, we find in terms of the quantities introduced above that ( J) 
in the direction of the nearest clump is of the order 

where the density profile inside the clump was considered roughly constant and the angular 
acceptance of the detector AQ was supposed to be greater or equal to the field of view of 
the clump AVld ~ {f/5y^^. We consider as an example the same choice of parameters as 
in Ref. 0: / ~ 0.01, M^i ~ IO^Mq, 6 ~ lOl In this case we find AVl^i ~ 4 ■ 10"^ sr, and 



taking AQ = 10^^ sr, we obtain ( J (ipd) ) ~ 3-10^ which we can compare with the analogous 
quantity in a smooth halo scenario. In Fig. 7 of Ref. the values of ( J ) {AQ = 10~^ sr) 
for a detector with AQ = 10~^ sr pointing towards the galactic centre are displayed (see 
also our Fig. |I|); the value of ( J (ipd) ) is about one half of the value for the most singular 
Navarro et al. profile, which gives detectable 7-ray lines for a relevant portion of MSSM 
models. 

The clump in our example might therefore be a very bright dark matter source, and 
a signal from neutralino annihilation into monochromatic photons in its direction could 
potentially be detected with an Air Cherenkov Telescope (ACT). In practice, the probability 
is small of detecting such a signal randomly pointing an instrument with a small angular 
acceptance. It might be of some help to combine ground- and satellite-based observations. 
A satellite detector, which has a wide field of view but also a much smaller effective area 
with respect to an ACT, may search in the continuum 7-ray spectrum for brighter spots 
in the sky which have no luminous counterpart. Such signals might then be identified as 
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clumps of dark matter if one would detect with an ACT the 7-ray lines from neutralino 
annihilations. For such a method to be practical, higher overdensities S may be needed. It 
should also be kept in mind that Eq. ( ]Tl| ) gives just a qualitative feature of the possible 
result; the possibility for the nearest clump of being much further away or a more realistic 
density profile may change that result by orders of magnitude. 

Much firmer predictions may be formulated in the many small clumps scenario; in this 
case we assume that most of the clumps cannot be resolved even by a detector with a rather 
small angular acceptance, say about AQ ~ 10~^ sr. There might still be some clumps which 
are resolvable just because they happen to be nearby and these should be treated as in the 
previous case. 

From Eq. (||), the probability for a clump of being at a line of sight distance (L, L + dL), 
a viewing angle defined by {cosip, cosip + dcosi/j) and at some azimuthal angle with respect 
to the direction of the galactic center (0, + dcp) is given by 

Pci{L,ip) dLdcosipdcp = —— p{L,ip) L? dLdcosipdcp . (12) 

Assuming that the clumps can be regarded as point-like sources, we can derive the analogue 
of Eq. (p!OD (as in the latter we factorize out 1/Att): 

} ]V r r 
{J{ij)){AQ) = ——-^^ / dQ' / dLp,i{L,ij')- 



8.5 kpC AQ JaQ J line of sight 



LW^^^' U.3GeV/cm3j • (13) 



Taking Eqs. (|^) and (|T2D into account, this can be rewritten as 

I fS f Po \ 
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8.5 kpc AQ V0-3GeV/cm 

dn' I dL ( f^"^''^ ] . (14) 

An J line of sight \0.3 GeV/cm'^ / 



Comparing Eq. (|1^ with Eq. ([T0|), we realize that in a scenario with many small clumps 
the angular dependence of the signal is different from the one in the smooth halo sce- 
nario. This is illustrated in Fig. |T| for two different halo models, picking as an example 
f6 = 20: a Navarro et al. profile (Eq. (|l]) with (a,/?, 7) = (1,2,1) and, in our exam- 
ple, Po = 0.3 GeV/cm^, a = 9 kpc) and a modified isothermal sphere (Eq. (|1]) with 
(a,/3,7) = (2,2,0), Po = 0.3 GeV/cm^, a = 3.5 kpc). The parameter fS mainly de- 
termines the relative importance of the smooth and clumpy components. An interesting 
feature, shown in the Figure for the Navarro et al. profile, is a possible break in the angular 
spectrum. This could be a possible signature to discriminate the signal from neutralino 
annihilations into continuum 7-rays from the galactic 7-ray background, and may be indeed 
suggested by a recent analysis of EGRET data . 

We are now ready to give predictions for the 7-ray flux from neutralino annihilations. To 
minimize the impact of the halo model and of experimental uncertainties, we concentrate 
on the flux at high latitudes, b > 60° and 0° < ^ < 360° {AQ = 0.84 sr), rather than 
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FIG. 1. The value of ( J{ip) )) (Afi) for two different halo profiles. The contribution from the 
smooth and clumpy component are also given. 



considering the flux towards the galactic centre which as shown in Fig. ^ is maximal. The 
modified isothermal profile of Fig. |I| gives 

( J(90°) y^""'^ (0.84sr) + ( J(90°) )"'™ (0.84sr) ^ 0.93 ■ (1 + 1.8 • f5) . (15) 

For simplicity we have made the reasonable assumption that / is small. If that is not true 
we have to replace 1 by (1 — /)^ in the above equation (as well as Eq. ( POj ) below). The 
analogous estimates with any of the halo models considered in Ref. are within a factor 
of 2 of the value given in Eq. (|T5|) . There is therefore a very weak halo model dependence 
in these results. In Fig. ^ (a) we plot the integrated 7-ray flux above the energy threshold 
Eth = 1 GeV for our set of MSSM models in the smooth halo scenario. Also shown in the 
figure is the corresponding 7-ray flux measured by the Energetic Gamma Ray Experiment 



Telescope (EGRET) as inferred from the analysis in Ref. [^ : 

^-^{E > 1 GeV) = (1.0 ± 0.2) x 10"*^ cm^^ g-^ sr"^ . (16) 

We can compare with this value to obtain a constraint on the allowed values of the parameter 
fS. It is however useful to analyse this together with the analogous constraint we can 
derive in the scenario of many small clumps from neutralino annihilations into cosmic ray 
antiprotons. 



B. Antiprotons 

Neutralino annihilations of relic neutralinos in the galaxy may produce cosmic ray an- 



tiprotons ( |1[] and references therein, p2| , p8| ) mainly from jets, in a process which is anal- 
ogous to the case of continuum 7-rays. To model the fragmentation process and extract 
information on the number and energy spectrum of the antiprotons produced we have again 



8 



10 



L Bergslmm. J. Edsja. P. Gondola and P. Ullio, 15 

Smooth halo 



EGRET measured y flux 



Gaugino-like 
Mixed 

Higgsino-like 



E*; = 1 GeV _ 



giO XXOXOOOB! 



ooia]ao+ 

O [Si's! 

I- I I 



10 



10 



10 



^ 10 
■> 

0) -6 

10 

} -7 

■<o 10 

'£ -8 

S 10 

■ Q. 

e 

-9 

10 



10 



-10 



10 



-11 



10 



-12 



10 



I. Bergslmm. J. Edsja. P. Gondola and P. Ullio, 1 



XXX 



10 



BESS 95 



Gaugino-like 
lUlixed 

Higgsino-liite 



Smootii iialo _ 
T= = 0.4 GeV 



10 10 
Neutralino Mass (GeV) 



Neutralino Mass (GeV) 

FIG. 2. The signal of (a) continuum gamma and (b) antiprotons versus the neutrahno mass. 
Only models with 0.025 < ^x^"^ < ^-^ have been included in this and the following figures. 



used the Lund Monte Carlo Pythia 6.115 and applied the same tabulation technique as 
for the production of photons. Including the same set of final states and treating the Higgs 
bosons in the same way, for any given MSSM model we can then obtain the energy spectrum 
of antiproton dark matter sources. 

If we assume a smooth distribution of WIMPs in the galaxy, the production rate of 
antiprotons in the volume element d^x at the galactic position x is given by 

— — d X = \ —— \ ■> vcTp ,^ d X (17) 

dT \M^j y ' dT ^ ' 

where T is the kinetic energy of the antiprotons. We will not discuss here the few-clumps 
scenario as those predictions are extremely sensitive to the parameters which define the 
model. We focus instead on the many small clumps scenario, treating again single clumps 
as point-like sources. In this case we find: 



dT 



d^x = NciPdix) J d^rcipliini) - ^vap d^x 

F 



It is not straightforward to simulate the propagation of charged particles in the galaxy. 
Different models have been proposed, and no consensus has been established yet. We present 
results in the limit in which the propagation is modeled by pure diffusion, using the analytic 



solution derived in Ref. [^] to which we refer for further details. 

The BESS collaboration has in a recent measurement of cosmic ray antiprotons 
found that the spectrum for the antiproton flux versus the kinetic energy T is consistent 
with being flat for T in the range between 180 MeV and 1.4 GeV. We consider the value of 
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the measured flux at some low value of the kinetic energy, where the 'trivial' antiproton flux 
generated by cosmic-ray reactions in the interstellar medium is believed to be less dominant. 
At T = 400 MeV the result found by BESS is 

^p{T = 400 MeV) = lAtU x 10"^ p cm^^ sr^^ GeV"^ . (19) 

In Fig. (b) we compare this value with the predictions for antiprotons from neutralino 



annihilations in a smooth halo scenario (i.e. the source given as in Eq. (17)) at the same 
energy, using for the diffusion model the same set of parameters as in Ref. [^], consider- 
ing appropriate values of the solar modulation parameters and picking as halo profile the 
modified isothermal distribution. It is indeed tempting to conclude that some of our models 
are already excluded by the BESS measurement. However, one has to keep in mind the big 
uncertainties involved, mainly in the antiproton propagation; for instance it is not clear how 
large a fraction of antiprotons generated in the halo can penetrate the wind of cosmic rays 



leaving the disk ||29|. We introduce in the flux predictions a rescaling factor k which contains 
the uncertainties deriving from the choice of the parameters which define the propagation 
model considered and from possible deviations from this simple approach. 

We consider now the many small clumps scenario. The production rate of antiprotons 
in this case is given by Eq. (0); the strength of the signal compared to the smooth case is 
again mainly determined by the product fS. At T = 400 MeV and for the same halo profile 
considered above, we find: 

$p- = k{l + 0.75 ■ f6) ■ (20) 

We have checked that the coefficient 0.75 depends very weakly on the halo profile considered 
and on T. A conservative limit on the dumpiness parameter f6 can be obtained choosing 
the uncertainty factor k as 

A; G [0.2,5]. (21) 



10 



^ 10 
1> 

<U -6 

^ 10 

-7 

■<o 10 

CM 

E .8 

3 10 

■ a 

-9 

10 



L. Bergstrom. J. Edsjb. P. Gondoio and P. Ultia, I 



T-p = 0.4 GeV 
^ E"; = 1 GeV 



10 



-10 



10 



-11 



10 



-12 



BESS 95 



Gaugino-like 
Mixed 

Higgsino-like 



Smooth halo 



-12 -11 -10 -9 -8 -7 -6 -5 

10 10 10 10 10 10 10 10 

«'cont.y(cm"'s-' sr-^) 



FIG. 4. The antiproton flux versus the continuum 7 flux for a smooth halo. 

We consider a value of f6 excluded if the whole range of possible antiproton fluxes given by 
Eqs. ([20| ) and (^11) exceeds the measured value, Eq. (|l^). 



C. Determining the dumpiness factor f5 

We have shown that in the many small clumps scenario the signals from dark matter 
annihilations into 7-rays and antiprotons depend critically on the dumpiness parameter fS. 
Focusing on the MSSM, we use the rescalings derived in Eqs. (15) and (12( 



to determine for 

each super symmetric model the maximal value of f6 for which the experimental constraints 
on the fluxes of continuum photons and antiprotons are not violated. This is shown in 
Fig. 1^ (a), where the maximal rescaling is given versus the neutralino mass, and where we 
use different symbols to indicate which of the two bounds is more restrictive. As can be seen, 
the antiproton flux puts the highest constraints on the dumpiness at low masses, whereas 
the continuum gammas put better constraints at higher masses. We see that the present 
experimental limits constrain fS < 10^ for all masses. 

As shown in Fig. ^, the two signals are strongly correlated since they are both produced 
from jets. In this sense the information we get from the two experimental limits is not entirely 
complementary. At higher masses, both fluxes go down since they are both proportional 
to correlation also decreases since the antiproton fluxes are only given in a 

small energy interval while the gamma ray fluxes are integrated above a threshold. Hence 
the antiproton flux in a given low energy interval decreases more than the gamma ray flux 
as we go to higher neutralino masses. In Fig. ^ (b) we analyse how restrictive one detection 
method is compared to the other. 

Having derived for each of the MSSM models in our sample the maximal allowed value 
for the dumpiness parameter f6, in the next section we analyse the consequences of this 
result for other indirect detection methods of neutralino dark matter. 
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FIG. 5. The number of expected photons in 0.84 sr towards b = 90° cohected in 4 years from 
(a) the 27 and (b) the Z7 final states. The expected 5a hmit from the GLAST detector is also 
shown. 



V. OTHER DETECTION METHODS 

In this section we consider the many small clumps scenario with the highest possible 
value of f6 as given in the previous section and investigate what effect that has on other 
dark matter searches. We fix again as smooth halo distribution to compare with the modified 
isothermal distribution, Eq. (0) with (a,/3,7) = (2,2,0), po = 0.3 GeV cm~^, a = 3.5 kpc 
and -Ro = 8-5 kpc. 



A. Monochromatic 7-ray lines 

As we have seen, the same scaling applies to the continuum and the line 7-ray signals, 
it is therefore straightforward to derive the maximal fluxes of monochromatic photons from 
neutralino annihilations. 

We perform this analysis in light of the potential of the next generation of satellite-based 
7-ray detectors, and in particular of the proposed Gamma-ray Large Area Space Telescope 
(GLAST) |3T|. To prevent uncertainties due to the choice of the dark matter halo profile to 
play any role in the following discussion, we fix again as field of view a 0.84 sr cone in the 
direction b = 90°. In the actual experiment the detector will collect data with a Air sr angular 
acceptance; as for most halo profiles the ratio signal to squared root of the background is 
greatly enhanced towards the galactic centre, the predictions we show are an underestimate 
of the possible results. 

Taking into account the screening of the earth, the useful geometrical acceptance of 
GLAST towards a fixed point of the sky in a 0.84 sr cone is 0.21 m^ sr the energy 



resolution is assumed to be 1.5%. We display in Fig. ^ the number of expected 7s in 4 years 
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FIG. 6. The difference of the diffuse neutrino flux towards the galactic centre to that to the 
antigalactic centre. The fluxes are averaged over 2.5 sr which maximizes signal to noise and they 
are rescaled maximally as allowed by the antiproton and continuum 7 fluxes. The limits are for a 
neutrino telescope with an exposure of 10 km^ yr. 



of exposure time when the fluxes have been maximally rescaled according to Fig. Also 
shown is the curve giving the minimum number of events needed to observe an effect at the 
5(7 level, where, in lack of data, we have assumed above 1 GeV a 2.7 power law falloff for the 
diffuse 7-ray background and inferred its normalization from Ref. p^]. As can be seen, a 



fair fraction of our set of supersymmetric models can be probed under these circumstances. 
Remember that the number of photons given in Fig. |^ is towards b = 90° and, depending 
on halo proflle, we expect more events towards the galactic centre, with a larger portion of 
the MSSM parameter space which might be probed. 



B. Diffuse neutrinos 

A possibility to get a detectable neutrino flux from WIMP annihilations which has rarely 
been considered in the literature is neutrinos from annihilation in the galactic halo. Of 
particular importance is XX ~^ W~^W~, since the W bosons decay in 10 % of the cases 
directly to a muon plus a muon neutrino with a hard neutrino spectrum, which may facilitate 
detection in neutrino telescopes. 

This flux would scale in exactly the same way as the gamma flux in the presence of 
clumps and with future O (km^) neutrino-telescopes, the diffuse neutrinos might prove 
more constraining than antiprotons and continuum 7s at high masses (several hundred GeV 
- TeV region) where the rescaling can be high {f6 > 10^). 

The flux has been calculated in essentially the same way as for neutralino annihilation 
in the Sun/Earth with the help of the Lund Monte Carlo Pythia 6.115. The only 
difference is that some annihilation products will decay and produce neutrinos in the halo. 
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whereas they are stopped before they decay in the Sun/Earth. 

The neutrino-induced muon flux from neutrahno annihilations in a smooth halo is about 
10~^-1 km~^ yr~^ sr~^ above 100 GeV. Compare this with the atmospheric background of 
about 9500 km~^ yr^^ sr~^ vertically and 30000 km^^ yr^^ sr^^ horizontally for this 



threshold. To be able to distinguish the signal from the background we have to rescale the 
fluxes by the allowed dumpiness factor derived in the previous section and we also have to 
make use of the fact that the signal is enhanced towards the centre of the galaxy. 

The best prospects are probably given by large-area neutrino telescopes with relatively 
high detection thresholds. We can imagine measuring the flux in a solid angle AQ towards 
the galactic centre and compare with the flux in the same solid angle in the opposite direction. 
The limit we can put on the flux is at the 3cr-level approximately given by 



[$,(0°)-$,(180°)],^,^3^^. (22) 

where S is the exposure. For the modified isothermal sphere, it turns out the best limits are 
obtained with = 2.5 sr, for which we obtain (J(0°))(Afi) = 4.16 and ( J(180°))(Afi) = 
1.09. In Fig. ^ we show the difference of the diffuse neutrino flux towards the galactic centre 
to that in the opposite direction for a muon energy threshold of 100 GeV. Also shown are 
the limits that can be reached with an exposure of 10 km^ yr. For different exposures, the 
limits change as the square root of the exposure. If we increase the threshold from 100 GeV, 
we can gain a small factor in sensitivity at higher masses, but lose at intermediate masses. 

An ideal neutrino detector for this signal would view the galactic center through the 
center of the Earth (i.e. it should be at 29 degrees latitude), since then the atmospheric 
background is minimal. The strength of the signal of course depends on the halo profile, 
but it is more likely that the halo profile is steeper towards the galactic centre than the 
isothermal sphere and hence the signal is even bigger then envisioned here. We might have 
to worry about other sources of high-energy neutrinos at the galactic centre (like neutrinos 
from the black hole believed to exist in the centre). These other sources can probably be 
removed by not looking at the very centre of the galaxy. 



C. Positrons 



From neutralino annihilation in the halo we would also get a flux of positrons which 
might be detected by satellite |Q or high-flying balloon experiments |^ 



The propagation 

of positrons is a more difficult issue than for antiprotons since positrons are so easily deflected 
and destroyed. We have calculated the positron fluxes using Pythia 6.115 and have 
used the propagation model in Ref. ||3^ with an energy dependent escape time (a more 



detailed investigation is in preparation |H^). In Fig. [7] we show the positron fluxes versus 
the neutralino mass when they have been rescaled with the maximal f6 allowed by the 
antiproton and the continuum gamma fluxes. We compare with the measurement by the 
HEAT experiment at 10 GeV. Also shown is the prediction of the background at this 
energy as given in Ref. [^. It would seem that the positrons put more stringent bounds 



on fS than the antiprotons and continuum 7s at higher masses. The positron fluxes are 
however even more uncertain than the antiproton fluxes and can easily be wrong by a factor 
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FIG. 7. The positron fluxes rescaled maximally as allowed by the antiproton and continuum 7 



fluxes. The 1994 HEAT measurement at 10 GeV [35| is shown together with the Strong et al. |3S] 
prediction for the background at this energy. 



of 10. Hence we can't use the positrons to constrain f5 further, but we see that we might 
be able to get measurable fluxes. 



D. Direct detection 

For the direct rates, we have used the procedures in Ref. Since these rates only 



depend on the local halo density at present, they will as expected not put any severe con- 
straints on the dumpiness of the halo as a whole. They will of course be much enhanced if we 
happen to be inside a clump at present. As with the neutrinos from neutralino annihilation 
in the Sun/Earth we can however have correlations between these signals and the signals 
giving high p or 7 fluxes. These correlations are not very strong, however. 



E. Neutralino annihilation in the Earth/Sun 

Neutrinos are produced in the annihilations through the decays of quarks, leptons and 
gauge bosons produced in the primary annihilation process. During the several billion 
years the Earth and Sun have existed, there may have been a substantial accumulation of 
neutralinos due to capture, i.e. scattering and subsequent gravitational trapping. 

The fluxes of neutrino- induced muons from neutralino annihilation in the Earth/Sun are 
mostly determined by the capture rates, which in turn depend on the local halo density. 
They are thus insensitive to different halo proflles; if the halo is clumped, there can however 
be fluctuations in the capture rates by time, but on the average we will capture the same 
amount of neutralinos as without clumps. The amount of these fluctuations in capture rate 
and consequently in annihilation rate depends on the time between encounters, the size of 
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the clumps, r,.|, and their overdensity 5. For the small-clumps scenario these fluctuations 
are expected to be small. 

Since these rates do not depend strongly on the clumping, they will not put better 
constraints on the dumpiness than the antiproton fluxes or the continuum gammas. 

VI. CONCLUSIONS 

To conclude, we have found that the limits on antiproton and continuum 7 fluxes already 
constrain in a non-trivial way the dumpiness of the Milky Way dark matter halo (if made of 
neutralinos). The general pattern is that at lower masses the p flux puts the best limits, but 
at higher masses the continuum 7 flux is better. Within the MSSM, the allowed dumpiness 
is less than fS ~ 10^ for all neutralino masses (assuming that the neutralinos make up most 
of the dark matter of our galaxy). 

We have also investigated what the detection prospects would be for other dark matter 
searches in this clumpy scenario, where the maximal rescaling is given by the limits on the 
antiproton and continuum 7 fluxes. We have found that the fluxes of monochromatic 7- 
lines from halo annihilation into the Zj and 27 final states can be enhanced enough to be 
seen by upcoming experiments like GLAST. We have also found that the fiux of positrons 
from neutralino annihilation in the halo gets high enough to even exceed the current limits 
from the HEAT experiment. The uncertainties for the positron flux arc particularly large, 
however, and at present we can merely conclude that it is possible to obtain measurable 
fluxes of positrons. The rarely-mentioned diffuse neutrino flux from neutralino annihilation 
in the halo can, for heavy neutralinos, be boosted enough to show a detectable difference in 
flux towards the galactic centre and the galactic anticentre. 

It is reassuring that new detectors, like GLAST for gamma rays and AMS for antiprotons, 
will obtain more stringent bounds on the dumpiness of the Milky Way halo. And, of course, 
finding evidence for neutralino annihilations in the halo would be one of the most important 
scientific discoveries of our time. 
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